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TECHNICAL  NOTE  3511 


DESCRIPTION  AND  ANALYSIS  OF  A  ROCKET- VEHICLE 
EXPERIMENT  ON  FLUTTER  INVOLVING  WING 
DEFORMATION  AND  BODY  MOTIONS1 
By  H.  J.  Cunningham  and  R.  R.  Lundstrom 


SUMMARY 


Flight  tests  and  a  mathematical  analysis  were  made  to  demonstrate 
and  confirm  a  type  of  subsonic  flutter  involving  rigid-body  motions  and 
wing  deformations.  For  the  configuration  considered,  the  period  of  the 
oscillation  was  approximately  100  chords  per  cycle  which  is  well  within 
the  range  of  period  found  in  dynamic-stability  work  on  rigid  aircraft 
with  free  controls.  A  mathematical  analysis  based  on  two-dimensional 
incompressible  flow  provided  a  conservative  prediction  of  the  airspeed 
at  which  the  low- frequency  flutter  occurred.  It  was  found  that  wing 
bending  stiffness  is  the  important  parameter  for  preventing  such  flutter. 


INTRODUCTION 


Interaction  of  deformations  of  an  aircraft  structure  with  the  passing 
airstream  can  lead  to  the  dynamic  instability  known  as  flutter.  For 
bending-torsion  wing  flutter,  the  frequency  of  oscillation  is  fairly  high 
and  usually  approaches  the  natural  torsional  frequency  of  the  wing  in 
still  air.  Such  an  oscillation  may  be  contrasted  with  ordinary  dynamic- 
stability  phenomena  involving  rigid-body  modes  which  lead  to  much  lower 
frequencies  that  are  usually  controllable. 

The  fact  that  the  calculated  flutter  speed  may  be  modified  by  the 
addition  of  free-body  modes  has  been  recognized  for  many  years.  For 
example,  about  20  years  ago  it  was  found  analytically  (ref.  l)  that  body 
mobility  had  a  slight  favorable  effect  on  the  calculated  flutter  speed 
of  a  particular  configuration  typical  of  that  day.  The  problem  has  from 
time  to  time  been  reconsidered  in  both  American  and  British  literature 


-^Supersedes  the  recently  declassified  NACA  RM  1*50129,  "description 
and  Analysis  of  a  Rocket-Vehicle  Experiment  on  Flutter  Involving  Wing 
Deformation  and  Body  Motions11  by  H.  J.  Cunningham  and  R.  R.  Lundstrom, 
1950. 
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and  the  necessity  for  determining  any  potential  detrimental  effect  of 
special  configurations,  sweptback  wings,  higher  speeds,  and  higher  alti¬ 
tudes  has  lately  become  more  insistent.  A  recent  paper  by  Broadbent 
(ref.  2)  discusses  the  necessity  for  including  free-body  modes  in  the 
study  of  sweptback  wings. 

Controlled  experimentation  involving  free-body  modes  is  highly  desir¬ 
able  although  difficult  even  at  low  speeds  in  a  wind  tunnel .  In  refer¬ 
ence  3  Lambourne  gives  some  experimental  results  and  describes  various 
difficulties  encountered,  principally  the  difficulty  of  supporting  the 
model  so  that  actual  flight  behavior  is  sufficiently  well  simulated. 

The  National  Advisory  Committee  for  Aeronautics  has  been  engaged  for 
several  years  in  making  flutter  experiments  in  the  transonic  range  by 
the  use  of  modern  telemetering  techniques  with  rocket-powered  vehicles. 

Some  of  the  results  of  these  experiments  are  given  in  several  MCA  papers 
including  references  4  and  5*  Such  flutter  research  has  been  concerned 
almost  entirely  with  wing  bending- torsion  flutter.  Recently,  however,  in 
connection  with  some  of  these  experiments,  flutter  failures  were  obtained 
which  definitely  involved  much  lower  frequencies  than  those  obtained  in 
ordinary  bending-torsion  flutter.  A  preliminary  account  of  an  unexpected 
low-frequency  failure  is  given  in  reference  4,  and  two  other  low-frequency 
failures  have  been  observed  by  W.  T.  Lauten,  Jr.,  and  J.  M.  Teitelbaum 
during  rocket-propelled  model  tests  at  the  Langley  Pilotless  Aircraft 
Research  Station  at  Wallops  Island,  Va. 

It  was  decided  to  repeat  the  experiment  of  the  D  model  of  reference  4 
with  more  extensive  instrumentation  in  order  to  obtain  information  specific 
to  this  type  of  flutter.  The  present  paper  includes  a  description  of  and 
results  from  model  D  of  reference  It-  as  well  as  similar  and  more  comprehen¬ 
sive  material  for  the  repeat  experiment  on  the  model  which  has  been  desig¬ 
nated  model  E.  Results  of  an  analysis  developed  on  the  baBis  of  certain 
simplifying  assumptions  are  given,  and  a  comparison  is  made  between  analyt¬ 
ical  and  experimental  results . 


SYMBOLS 


Py 

a  nondimensional  position  of  wing  elastic  axis,  — —  -  1 

100 

a  +  %  nondimensional  position  of  wing- section  center  of  gravity, 


100 

Ag  geometric  aspect  ratio  of  one  exposed  wing  panel,  for  rectangular 

plan- form  wings,  Length/Chord 
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k 

l 
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M 

Mcr 


wing  semichord,  ft 
bending  rigidity,  lb- in. 2 

nondimens ional  distance  (in  wing  semi chords)  from  wing  midchord 
line  to  missile  center  of  gravity  measured  parallel  to  missile 
longitudinal  axis,  positive  when  center  of  gravity  is  behind 
wing  midehord  line 

first-bending  natural  frequency,  cps 
second-bending  natural  frequency,  cps 


first-torsion  natural  frequency,  cps 

uncoupled  first-torsion  frequency  relative  to  elastic  axis 


ft 


i  - 


(%/r4 


-  (Wft) 


1/2 


cps 


y 


structural  damping  coefficient  in  bending 

structural  damping  coefficient  in  torsion 

mathematical  quantity  having  qualities  of  structural  damping 
coefficient 

torsional  rigidity,  lb- in. 2 

moment  of  inertia  of  missile  about  center  of  gravity,  slug- ft2 

reduced  frequency  parameter,  aib/v;  it/ k  is  period  of  oscillation 
in  chords  per  cycle 

length  of  wing  along  leading  edge  outboard  of  body,  in. 

mass  of  wing  per  unit  length,  slugs/ft 

Mach  number 

theoretical  Mach  number  at  which  sonic  velocity  is  first  attained 
over  section  of  wing  taken  perpendicular  to  leading  edge  at 
zero  lift 


v 
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%  nondimens ional  radius  of  gyration  of  wing  section  about  elastic 

axis  faz/n ib2  where  is  the  mass  moment  of  inertia  of 

wing  about  its  elastic  axis  per  unit  length  in  slug-ft^/ft 


t 

v 

V 

btuh 

xcg 

xea 

K 

k ( stnd . ) 
P 


time  after  rocket  launching,  sec 
speed,  ft/sec 

nondimensional  flutter-speed  coefficient 

distance  of  center  of  gravity  of  wing  section  behind  leading 
edge,  percent  chord 

distance  of  wing  elastic  axis  behind  leading  edge,  percent 
chord 

mass  ratio,  npb^/m 

value  of  k  when  p  is  standard  sea- level  density 
air  density,  slugs/cu  ft 


Cd 

angular 

coh 

angular 

% 

angular 

2rtfa., 

fi 

flutter 

flutter  parameter,  (l  +  ig) 


EXPERIMENTAL  INVESTIGATION 
Apparatus  and  Methods 


The  test  vehicle  was  essentially  a  tailless  configuration  having  the 
test  wings  as  the  only  stabilizing  surfaces  in  pitch.  The  model  was 
powered  by  a  modified  Aerojet  12AS-1000  D  rocket  motor  capable  of  carrying 
it  to  greater-than-sonic  speed  with  an  acceleration  of  about  k  times  grav¬ 
ity.  A  sketch  of  model  E  is  shown  in  figure  1.  Center  of  gravity,  weight, 
and  moment  of  inertia  in  pitch  varied  with  flight  time  because  of  the  fuel 
consumed.  These  parameters  plotted  against  flight  time  for  model  E  are 
shown  in  figure  2.  The  launching  was  made  at  the  Langley  Pilotless  Aircraft 
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Research  Station  at  Wallops  Island,  Va.,  from  a  near- zero- length  launcher, 
as  described  in  reference  5-  Photographs  of  model  E  in  the  launching 
position  are  shown  in  figure  5.  The  test  wings  were  made  of  laminated 
white  pine  with  an  inlay  of  0.032-inch  24S-T  aluminum  alloy  to  duplicate 
as  closely  as  possible  the  test  wings  of  model  D  of  reference  4.  Strain 
gages  were  mounted  on  the  elastic  axes  of  the  wings  4  inches  from  the  wing 
root,  as  shown  in  figure  1. 

A  six- channel  telemeter  was  installed  in  the  model  with  instrumen¬ 
tation  to  give  continuous  readings  of  the  strains,  left-wing  bending, 
right-wing  bending,  and  right-wing  torsion,  at  the  corresponding  gage 
locations  and  also  to  give  continuous  readings  of  longitudinal  and  normal 
accelerations  of  the  model  center  of  gravity  and  angle  of  attack  of  the 
missile.  Speed  of  the  model  was  obtained  by  integration  of  longitudinal 
acceleration,  and  the  altitude  was  obtained  from  a  pulse-type  tracking 
radar  unit.  Atmospheric  conditions  prevailing  at  the  time  of  flight  were 
obtained  from  a  radiosonde. 

The  accuracy  of  the  quantities  measured  is  believed  to  be  within  the 
following  limits:  longitudinal  and  normal  accelerations,  2  percent;  vel¬ 
ocity,  4  percent;  wing  bending  and  twisting  moments,  20  percent;  and  phase 
angle  between  any  two  quantities,  15°. 


Experimental  Results 

The  physical  characteristics  of  the  wings  as  determined  by  preflight 
ground  tests  are  listed  in  table  I.  Data  for  model  E  are  presented  in 
figures  4  to  7  as  functions  of  the  flight  time.  Included  in  these  figures 
are  the  available  data  for  model  D  from  reference  4  for  comparison.  Fig¬ 
ure  4  includes  flight  velocities;  figure  5  shows  Mach  number  and  air  density; 
figure  6  shows  longitudinal  acceleration.  Figure  7  presents  the  following 
data:  bending  moments  of  both  wings  (positive,  wing  bent  down)  at  the 
strain-gage  locations,  twisting  moment  of  the  right  wing  about  its  elastic 
axis  (positive,  leading  edge  up)  at  the  strain-gage  location,  and  normal 
acceleration  of  the  missile  center  of  gravity  (positive,  up) .  It  was 
desired  to  obtain  the  angle  of  attack  of  model  E  directly  as  a  function 
of  time,  but  the  angle- of-attack  indicator  was  inoperative  during  flight. 

As  a  result,  the  angle  of  attack  was  determined  by  use  of  the  normal  accel¬ 
eration  and  an  assumed  lift-curve  slope  of  0.08  per  degree  as  obtained 
from  lifting- surface  theory. 

Comparison  of  figure  7  for  model  E  and  figure  9  of  reference  4  for 
model  D  indicates  that  both  flights  had  similar  behavior.  At  a  speed  of 
about  400  feet  per  second,  there  began  pitching  oscillations  with  slowly 
increasing  amplitudes  which  continued  to  a  speed  of  about  620  feet  per 
second  for  model  D  and  500  feet  per  second  for  model  E.  Over  this  interval 
the  primary  frequencies  of  oscillation  increased  from  4  to  7  cycles  per 
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second  for  model  E  and  from  5  to  8  cycles  per  second  for  model  D,  and 
there  was  a  considerable  component  of  the  wing  first-bending  natural  fre¬ 
quencies  present  on  the  wing  bending  strain-gage  traces.  An  interruption 
or  damping  of  the  oscillations  then  occurred  and  was  followed  by  final 
divergent  oscillations  with  frequencies  of  about  9-0  cycles  per  second 
for  model  D  and  8.5  cycles  per  second  for  model  E  which  resulted  in  fail¬ 
ures  of  the  wing  after  about  7  or  8  cycles.  The  final  oscillations 
were  of  a  much  more  regular  character  and  approached  a  divergent  sine 
wave.  The  speeds  at  the  beginning  of  the  final  divergent  oscillations, 

670  feet  per  second  for  model  D  and  580  feet  per  second  for  model  E,  are 
considered  to  be  the  experimental  flutter  speeds  because  of  the  regularity 
of  the  subsequent  oscillations.  As  indicated  by  the  time  histories,  the 
periods  of  the  low-frequency  oscillations  for  both  models,  from  the  begin¬ 
ning  of  a  detectable  oscillation  until  wing  failure,  were  about  100  chords 
per  cycle. 

The  primary  interest  in  flutter  research  centers  on  the  speed  and 
frequency  of  oscillation.  A  secondary  interest  exists  in  the  deflection 
amplitudes  and  in  the  phase  relationships  of  the  various  motions;  the 
latter  are  usually  not  experimentally  determinable  with  good  accuracy. 

From  analysis  of  the  oscillograph  records  of  model  E  taken  during  flutter 
over  the  interval  t  =  5*^0  seconds  to  5*^5  seconds,  the  amplitudes  of 
wing  deformation  and  body  motions  relative  to  a  wing- tip -bending  amplitude 
of  1  inch  were  found  to  average  about  O.OO65  radian  for  wing-tip  torsion, 
0.020  radian  for  pitching,  and  0.21  inch  for  vertical  translation.  The 
wing-tip-bending  and  wing-tip-torsion  amplitudes  were  determined  through 
the  assumption  of  a  constant  lift  distribution  from  wing  tip  to  wing  tip, 
whereas  the  amplitude  of  vertical  translation  was  obtained  by  integrating 
the  normal  acceleration  twice.  The  body-pitching  amplitudes  were  determined 
by  using  the  normal  acceleration  and  an  assumed  lift- curve  slope  of 
0.08  per  degree.  The  overall  accuracy  of  the  experimentally  determined 
amplitudes  is  thought  to  be  within  ±50  percent;  thus,  these  values  are 
at  least  of  the  proper  order  of  magnitude. 

It  can  be  seen  from  figure  7  that  the  phase  relationship  is  as  fol¬ 
lows  within  the  limit  ±15°  of  experimental  accuracy:  At  an  instant  when 
the  wing  is  bent  up  a  maximum  amount,  it  is  also  twisted  (leading  edge 
up)  a  maximum  amount  and,  as  interpreted  from  the  normal  acceleration, 
the  missile  is  pitched  nose  up  and  translated  down  a  maximum  amount. 


AMLYTICAL  INVESTIGATION 
Bases  of  Analysis 


Only  an  outline  of  the  analytical  treatment  used  to  obtain  the  results 
is  given  here.  It  is  assumed  that  the  theory  of  linear  superposition  for 
small  disturbances  holds.  A  simplified  configuration  which  has  four  degrees 
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of  freedom  in  a  symmetric  oscillation  is  treated.  There  are  two  rigid- 
hody  degrees  of  freedom:  pitching  about  the  center  of  gravity,  and  ver¬ 
tical  translation;  and  two  wing-deformation  degrees  of  freedom  which  are 
approximated  by  the  first-uncoupled-bending  and  first-uncoupled- torsion 
modal  shapes  of  an  ideal  uniform  cantilever  beam.  For  convenience,  the 
aerodynamic  coefficients  used  for  a  wing  section  are  those  of  incompres¬ 
sible,  two-dimensional  potential  flow  and  yield  forces  which  are  propor¬ 
tional  to  the  displacement  and  motion  of  that  section  (strip  analysis). 
Spanwise  distributions  of  air  forces  due  to  symmetric  motions  of  a  rigid 
body  are  thus  constant  in  magnitude  and  phasing  from  wing  tip  to  wing  tip 
including  the  body  intercept.  In  view  of  the  Mach  number  of  the  experi¬ 
mental  flutter,  which  is  well  below  the  critical  Mach  number  Mcr  of  the 

airfoil  section,  aerodynamic  coefficients  for  incompressible  flow  are 
considered  satisfactory.  Coefficients  for  compressible  flow  could  be 
substituted,  however . 

The  configuration  has  no  other  horizontal  airfoil  surface  and  no 
other  aerodynamic  forces  are  assumed  on  the  body.  Each  half -wing  is 
uniform  and  the  two  half -wings,  right  and  left,  of  each  missile  are 
treated  as  being  identical.  The  effect  of  gravity  on  an  oscillation 
having  the  frequency  of  the  observed  oscillation  is  considered  negligible 
as  are  the  effects  of  rocket  thrust. 


Equations  of  Equilibrium 

In  order  to  obtain  the  theoretical  conditions  for  flutter,  Lagrange  * 
equation  is  employed,  as,  for  example,  in  reference  6,  to  derive  four 
equations  of  equilibrium.  With  the  use  of  the  specified  condition  of 
harmonic  motion,  the  four  equations  may  be  written  in  matrix  form  as 
follows : 


ail  a12  a15  al4 

xl 

^X1 

a21  a22  a25  &2k 

x2 

“S)1* 

a3l  a32  a33  a34 

x3 

ftpxj 

a4l  a42  ali-3  a44 

3_ 

OTxj_ j. 

2 

In  this  equation  the  quantity  n  =  (l  +  ig)  where  g  has  the  prop¬ 

erties  of  a  structural  damping  coefficient.  (Structural  damping  force 
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is  proportional  to  amplitude  and  in  phase  with  velocity.)  The  quantities 
X]_,  x2,  x^,  and  xi,.  represent  the  nondimens ional  amplitudes  of  oscil¬ 

lation  in  wing  tending,  wing  torsion,  body  pitching)  and  vertical  trans¬ 
lation,  respectively.  On  the  basis  of  the  foregoing  assumptions,  for 
the  configuration  treated. 


all  ~  ~  ~  Ach 


a12  =  1.3558(2p  -  Aca) 

a13  =  1 . 5662  ^~s  .a-+-  ^  -  Acpj 

=  1.5662a1]_ 

*21  =  0.6778(^  -  Aah) 


a  -  ra2  A 

a22  "  ~Z - Aaa 


a25  =  1.2732 

a^  =  1  •  87  84a2]_ 

where 


ra2  +  (-s  +  a)*, 


a 


31  =  0>33177  (~S  +■  *  -  **  - 

ra2  .  (-s  +  a)% 


a52  =  0.5395 

rp2  . 

a33  =  —  '  AdP 
a3b  =  “Adh 
s-ki  =  0.33177an 

a42  =  0.39793a12 


"  A  da 


1 


ap  43 


a,.,  =  -A 


•cp 


=  “  Ach 


ACh  =  Ach(k)  =  -1  -  f  +  *¥ 


Aca  -  ^ax  ( a) 
^ah  -  ^-ahC^^3*) 
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Aaa  =  Aaa<k>a)  =  "  J  "  a2  "  +  a)p  +  (J  "  + 

-&♦*$] 

Acp  =  Aca(k> s) 

Adh  =  Aah(k,s) 

Adp  =  Aaa(k, s) 

Aap(k,a,s)  =  Aa a  -  (s  -  a)Aa^ 

Ado/k;a>s)  =  Aaa  "  (s  -  a)A 

ca 


and 


Mass  moment  of  missile  In  pitch  about  Its  center  of  gravity 

(Span)jtpb^ 


1  _  Total  mass  of  missile 
KT  (Span)jtpb2 


The  quantity  a  is  the  nondimensional  distance  (in  wing  semichords) 
from  the  wing  midchord  line  to  the  wing  elastic  axis,  positive  when  the 
elastic  axis  is  behind  the  midchord.  The  quantity  s  is  the  nondimen¬ 
sional  distance  (in  wing  semichords)  from  the  wing  midchord  line  to  the 
center  of  gravity  of  the  missile,  positive  when  the  center  of  gravity  is 
behind  the  midchord.  The  quantities  F  and  G  are  the  real  and  imagi¬ 
nary  parts,  respectively,  of  the  complex  function  C  =  C(k)  =  F(k)  +  iG(k) 

developed  in  reference  7 .  The  reduced  frequency  parameter  k  equals  — . 

It  is  to  be  recognized  that  equation  (l)  represents  a  characteristic- 
value  problem  of  aeroelastic  harmonic  motion.  The  a^'s  depend  on  the 

reduced  frequency  k,  whereas  Q  contains  the  frequency  w.  Combinations 
of  ca  and  k  which  result  in  a  specified  value  of  g  are  the  character¬ 
istic  values  (eigenvalues)  of  the  system.  Knowledge  of  k  and  co  leads 
to  a  flutter  speed  v.  The  trivial  solution  xq  =  X2  =  X3  =  xk  =  0  is  of 
course,  not  sought.  J  ’ 
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The  quantities  p  and  T  are  proportional  to  any  mechanical 
restoring  force  due  to  tody  motions  in  the  tody-pitching  and  vertical- 
translation  degrees  of  freedom,  respectively.  Such  restoring  forces  could 
exist,  for  example,  if  the  model  were  mounted  on  springs  in  a  wind  tunnel. 
When  the  wing-body  combination  is  of  the  freely  flying  type,  however,  there 
is  no  mechanical  restoring  force;  that  is,  the  natural  (in  vacuo)  frequen¬ 
cies  of  the  two  rigid-tody  degrees  of  freedom  are  zero.  Hence,  the 
characteristic- value  solutions  for  equation  (l)  are  sought  at  the 
limit  as  p  and ■  T  approach  zero . 


DISCUSSION  AND  COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  RESULTS 
Rayleigh-Bit z  Treatment  in  Four  Degrees  of  Freedom 


When  equation  (l)  for  harmonic  motion  is  treated  in  the  Rayleigh- 
Ritz  manner  there  exist  in  general  four  roots  of  flutter  speed  and  fre¬ 
quency  which  satisfy  the  equations  of  equilibrium,  although  one  or  more 
roots  may  have  no  physical  significance.  For  the  configuration  tested 
and  analyzed,  the  critical  (lowest)  flutter  speed  root  (designated  root  A) 
corresponds  to  an  oscillation  involving  appreciable  proportions  of  the 
rigid-body  stability  modes  at  a  frequency  which  is  a  fraction  of  the  wing 
first-bending  natural  frequency.  The  next  higher  speed  root  (designated 
root  B)  corresponds  to  the  conventional  wing  flexure- tor si on  flutter  and 
only  minute  amplitudes  of  rigid-body  motions  are  present. 

During  the  flight  time  after  launching,  several  flutter  parameters 
were  continually  changing.  Such  varying  parameters  were:  air  density, 
weight  and  mass  moment  in  pitch  of  the  missile,  and  the  location  of  the 
missile  center  of  gravity.  These  combined  changes  have  only  a  moderate 
effect,  an  increase  of  about  7  or  8  percent,  on  the  calculated  flutter 
speed  (root  A)  from  t  =  0  to  the  time  of  wing  failure.  The  calculated 
flutter  speeds  as  functions  of  time  are  included  in  figure  4.  (The 
corresponding  changes  in  flutter  frequencies,  amplitude  ratios,  and 
phasings  are  also  small.)  The  intersections  of  calculated  flutter  speeds 
and  flight  speeds  give  the  predicted  flutter  speeds,  520  feet  per  second 
for  model  D  and  435  feet  per  second  for  model  E,  as  shown  in  figure  4. 

The  corresponding  predicted  flutter  frequencies  are  6.0  cycles  per  second 
for  model  D  and  4.3  cycles  per  second  for  model  E.  Each  configuration 
did  not  flutter,  however,  until  its  flight  speed  was  30  percent  higher, 

67O  feet  per  second  for  model  D  and  58O  feet  per  second  for  model  E,  as 
also  shown  in  figure  4.  Thus,  the  analysis  was  conservative  in  its  pre¬ 
dictions  of  flutter  speeds. 

With  regard  to  the  degree  of  conservatism  of  the  analytically  pre¬ 
dicted  flutter  speeds  it  should  be  pointed  out  that,  even  though  the 
aspect  ratio  is  7  (including  body  intercept),  the  finite-span  correction 
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to  the  flutter  speed  is  large  for  small  values  of  the  reduced  frequency 
parameter  k.  For  models  D  and  E,  the  analytical  value  of  k  was  approx¬ 
imately  0.027  which  is  small  for  flutter  work.  Reference  8  gives  some 
results  of  flutter  tests  made  in  part  to  determine  the  effect  of  finite 
aspect  ratio  on  the  flutter  of  cantilever  wings .  The  flutter  was  the 
tending-torsion  type  and  occurred  mostly  in  a  range  of  k  from  about 
0.15  to  0.30.  According  to  the  results  of  reference  8  for  an  aspect 
ratio  of  6  to  8,  the  finite  span  causes  an  increase  of  roughly  5  to  10  per¬ 
cent  in  the  flutter  speed.  For  such  a  low  value  of  k  as  was  obtained 
with  models  D  and  E,  it  is  known  that  the  finite  span  causes  a  consider¬ 
ably  larger  increase  in  the  flutter  speed.  Therefore,  it  appears  that, 
if  the  proper  finite-span  corrections  were  made,  an  analytically  predicted 
flutter  speed  would  agree  much  better  with  the  experimentally  determined 
one.  Compressibility  effects  also  influence  the  results  and  can  be  taken 
into  account.  It  is  of  interest  to  note  that  an  oscillation  with  k  =  0.027 
has  a  period  of  117  chords  per  cycle  since  jt/k  equals  the  period.  This 
value  of  the  period  is  well  within  the  range  of  period  treated  in  dynamic- 
stability  work  with  rigid  or  undeformable  aircraft. 

The  amplitude  ratios  and  phasing  of  the  various  degrees  of  freedom 
associated  with  root  A  for  model  E  when  t  =  5*2  seconds  based  on  a  wing- 
tip-bending  amplitude  of  1  inch  are  as  follows: 


Degree  of  freedom 

Relative  amplitude 

Relative  phase  angle, 
deg 

Wing-tip  bending 

1.0000  inch 

0 

Wing-tip  torsion 

.00^36  radian 

179.1 

Rigid-body  pitching 

.01679  radian 

186.1 

Rigid-body  vertical 
translation 

.2^37  inch 

186.0 

The  parameters  used  were:  mass  moment  of  missile  in  pitch,  13.25  slug- 
feet2;  weight  of  missile,  211  pounds;  air  density,  0.002260  slug 
per  cubic  foot;  and  the  parameters  of  the  right  wing  of  model  E.  The  wing— 
tip-torsion  amplitude  is  one-quarter  of  the  body-pitching  amplitude,  and 
the  vertical  translation  amplitude  is  one-quarter  of  the  wing- tip-bending 
amplitude.  Each  degree  of  freedom  is  very  nearly  in  phase  or  almost  dia- 
,  metrically  out  of  phase  with  other  degrees  of  freedom,  (on  the  basis  that 

positive  translation  and  wing  bending  are  downward,  whereas  positive  angu¬ 
lar  displacement  is  leading  edge  or  nose  up).  Also  of  interest  is  the 
*  virtual  identity  of  phasing  of  the  two  rigid-body  modes,  which  fact,  together 
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with  the  appropriate  amplitude  ratios,  indicates  that  the  missile  is 
pitching  effectively  about  a  lateral  axis  l4.9  inches  forward  of  the  mis¬ 
sile  center  of  gravity.  All  the  analytical  amplitude  ratios  and  phase 
relationships  agree  with  their  experimental  counterparts  (reported  in 
section  entitled  "Experimental  Results")  within  the  experimental  accuracy. 

The  amplitude  ratios  of  wing  flexure-torsion  flutter  (root  B),  based 
on  a  wing- tip-bending  amplitude  of  1  inch,  are  minute  for  the  two  body 
degrees  of  freedom,  whereas  the  wing-tip-torsion  amplitude  is  O.O676  radian 
with  a  phasing  of  about  -l8°.  The  ratio  of  body  weight  to  wing  weight 
ranged  from  38  at  time  of  launching  to  31  at  time  of  wing  failure  for 
model  E;  if  this  ratio  were  smaller,  more  relative  body  motion  might  be 
associated  with  root  B. 


Effect  of  Various  Binary  and  Ternary  Combinations  of  Freedoms 

Analyses  employing  various  binary  and  ternary  combinations  of  degrees 
of  freedom  were  made  and  the  results  were  compared  with  those  of  the  four- 
degree-  of-  freedom  analysis.  Results  are  given  in  table  II.  The  parameters 
used  throughout  are  those  of  the  right  wing  of  model  D  and  those  listed  at 
the  top  of  table  II.  The  four-degree- of- freedom  analysis  is  designated 
type  (a). 

Consider  root  A.  The  addition  of  the  wing-torsion  degree  of  freedom 
to  an  analysis  not  already  including  it  corresponds  to  a  reduction  of 
torsional  stiffness  from  an  infinite  value  to  the  finite  experimentally 
determined  value.  Even  though  the  amplitude  of  wing-tip  torsion  in  a 
four- degree- of- freedom  analysis,  type  (a),  is  small  compared  with  body¬ 
pitching  amplitude,  addition  of  wing  torsion  to  analysis  types  (c)  or  (e) 
so  that  they  become  type  (d)  or  (a),  respectively,  causes  a  significant 
reduction  of  about  13  percent  in  the  analytically  determined  flutter 
speed. .  These  reductions  indicate  that  it  is  potentially  dangerous  to 
exclude  the  wing  torsion  from  an  analysis  of  body- freedom  flutter.  The 
addition  of  the  vertical- translation  degree  of  freedom  to  analysis 
type  (c)  or  (d),  so  that  analysis  type  (e)  or  (a),  respectively,  is 
obtained,  effects  an  increase  of  50  percent  in  the  flutter  speed  and 
demonstrates  clearly  that  it  is  not  realistic  to  exclude  vertical  trans¬ 
lation  from  an  analysis  treating  body  freedoms,  at  least  for  the  present 
configuration.  Thus,  it  seems  clear  that  at  least  the  four  degrees  of 
freedom  employed  in  this  report  should  be  included  in  an  analysis  of  sym¬ 
metrical  body- freedom  flutter.  The  wing-bending  stiffness  is  of  predomi¬ 
nant  importance  to  the  flutter  speed  for  root  A  inasmuch  as  that  flutter 
speed  is  virtually  directly  proportional  to  the  natural  bending  frequency, 
at  least  in  the  range  of  low  ratio  of  bending  frequency  to  torsion 
frequency. 

The  flutter  speed  and  frequency  of  root  B,  primarily  wing  bending- 
torsion  flutter,  are  unchanged  regardless  of  which  type  of  analysis,  (a), 
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(t>),  or  (d) ,  is  employed,  and  for  root  B  the  flutter 
affected  by  the  torsional  stiffness  in  contrast  with 
of  bending  stiffness  on  the  flutter  speed  of  root  A. 


Effect  of  Relative  Positions  of  Wing  and  Missile 
Center  of  Gravity 

Another  variable  whose  effect  was  studied  was  the  distance  from  the 
missile  center  of  gravity  to  the  wing  leading  edge.  The  flutter  parameters 
used  were  again  those  of  the  right  wing  of  model  D  and  those  listed  in 
table  II  except  for  the  fixed  distance  between  the  wing  leading  edge  and 
missile  center  of  gravity.  The  analysis  was  type  (a);  that  is,  it  included 
all  four  degrees  of  freedom. 

Figure  8  shows  the  two  calculated  roots  A  and  B  of  flutter- speed 
coefficient  as  a  function  of  the  distance  of  the  wing  leading  edge  behind 
the  center  of  gravity  of  the  missile.  The  experimental  flutter-speed 
coefficient  for  model  D  is  included  in  the  figure  at  a  wing  position  of 
3.38  semichords,  even  though  the  experimental  parameters  differed  slightly 
or  moderately  from  those  at  the  top  of  table  II  on  which  the  analytical 
curves  are  based.  If  the  parameters  of  the  analysis  coincided  with  their 
experimental  counterparts,  the  calculated  curves  would  be  affected  by  only 
a  few  percent.  The  important  point  is  that  the  experimental  flutter  speed 
is  about  30  percent  higher  than  the  calculated  one.  It  is  pertinent  to 
recall  that  the  missile  has  no  horizontal  airfoil  surface  other  than  the 
wing  and  that  the  missile  would  become  statically  unstable  with  the  wing 
very  far  forward  near  the  center  of  gravity.  Such  a  rigid-body  static 
instability  is  not  indicated  in  the  figure.  Two  low-frequency  failures 
of  missiles  which  did  have  horizontal  tails  and  which  had  the  missile 
center  of  gravity  at  approximately  one-eighth  chord  have  been  observed 
during  rocket-propelled  model  tests  made  at  the  Langley  Pilotless  Aircraft 
Research  Station  at  Wallops  Island,  Va. 

If  analysis  types  (c)  or  (d)  (see  table  II)  which  do  not  include  the 
freedom  of  vertical  translation  were  employed,  it  would  be  predicted,  on 
the  basis  of  the  aerodynamic  forces  of  two-dimensional  incompressible 
potential  flow  developed  in  reference  7>  that  single-degree-of-freedom 
rigid-body  pitching  instability  would  be  encountered  over  the  approximate 
range  of  wing  position  shown  in  figure  8.  This  single-degree  pitching 
instability  was  investigated  analytically  by  Smilg  in  reference  9-  The 

quantity  lo/rtpb^  of  reference  9  (where  Iq,  represents  the  mass  moment 
in  pitch  of  the  missile  divided  by  the  span)  is  about  ten  to  eleven 
thousand  for  models  D  and  E.  The  predicted  borderline  flutter  speed 
(root  A) ,  if  no  vertical  translation  were  permitted,  would  be  zero  in  the 
critical  range  of  wing  location,  since  no  mechanical  restraint  of  pitching 
exists;  that  is,  the  in- vacuo  natural  rigid-body  pitching  frequency  is 


speed  is  predominantly 
the  predominant  effect 
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zero.  It  appears,  however,  that  such  an  unhappy  possibility  is  precluded, 
even  in  two-dimensional  incompressible  flow,  by  the  actual  existing  free¬ 
dom  of  vertical  translation  of  the  missile  in  flight,  at  least  for  the 
configuration  treated. 

The  two  low-frequency  failures  of  missiles  observed  during  the  rodket- 
propelled  model  tests  mentioned  previously  will  now  be  considered.  Each 
of  these  models  had  a  pair  of  unswept  rectangular  uniform  model  wings 
2k  inches  by  8  inches.  For  convenience,  the  models  are  designated  models  6 
and  7.  The  wings  of  model  6  were  designed  to  be  similar  to,  although 
slightly  smaller  than,  the  wings  of  models  D  and  E  of  the  present  report, 
and  the  wing  first-bending  and  first-torsion  frequencies  of  model  6  were 
practically  identical  to  those  of  model  D.  The  leading  edges  of  the 
wings  of  models  6  and  7  were  located  slightly  ahead  of  the  missile  center 
of  gravity,  -0.25  semichord  for  model  6  and  -0.38  semichord  for  model  7* 

The  complete  models  6  and  7  had  weights  and  mass  moments  in  pitch  approx¬ 
imately  one-third  as  great  as  those  of  models  D  and  E.  Each  model  included 
a  horizontal  tail  having  an  area  equal  to  about  20  percent  of  the  wing  area, 
and  a  tail  length  of  about  one-quarter  of  the  total  tip-to-tip  wing  span. 

Thus,  it  may  be  seen  that  models  6  and  7*  except  for  the  added  hori¬ 
zontal  tails,  were  dynamically  fairly  similar  to  a  hypothetical  model  D 
or  E  which  has  had  its  wings  moved  forward  near  to  the  missile  center  of 
gravity.  In  the  light  of  such  observation,  it  is  interesting  to  examine 
figure  8  for  a  small  negative  value  of  the  abscissa  or  wing  position. 

In  that  region  the  speed  of  the  flutter  involving  body  motions  is  pre¬ 
dicted  to  be  lower  than  the  speed  for  wing  bending-torsion  flutter  of  a 
modified  model  D  or  E. 

All  these  circumstances  lead  to  the  conjecture  that  the  horizontal- 
tail  volumes  of  models  6  and  7  were  sufficiently  large  to  prevent  a 
static  instability,  but  not  sufficient,  or  the  tails  were  improperly 
placed,  to  prevent  a  low-frequency  body- freedom  flutter. 


CONCLUSIONS 


The  following  conclusions  may  be  drawn  from  flight  tests  and  a 
mathematical  analysis  made  to  investigate  a  type  of  flutter  involving 
wing  deformation  and  body  motions  of  a  particular  configuration  and 
certain  modifications  of  that  configuration: 

1.  The  possibility  of  a  flutter  type  of  dynamic  instability 
involving  appreciable  proportions  of  rigid-body  motions  as  well  as  wing 
structural  deformations  has  been  predicted  analytically  and  confirmed 
experimentally.  The  flutter  was  entirely  different  from  conventional 
wing  bending-torsion  flutter  in  that  the  frequency  was  about  one-half 
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the  wing  first-hending  natural  frequency.  The  period  was  approximately 
100  chords  per  cycle  which  is  well  within  the  range  of  period  found  in 
rigid-body  control-free  stability  work. 

2.  A  flutter  analysis  based  on  two-dimensional  incompressible  flow 
provided  a  conservative  (by  30  percent)  prediction  of  the  airspeed  at 
which  the  flutter  would  occur  for  the  configuration  studied. 

3-  Wing-bending  stiffness  rather  than  wing-torsional  stiffness  has 
the  predominant  effect  on  speed  of  the  body- freedom  flutter. 

k.  At  least  four  degrees  of  freedom,  wing  bending,  wing  torsion, 
body  pitching,  and  body  vertical  translation,  should  be  included  in  an 
analysis  of  symmetrical  flutter  involving  wing  deformations  and  body 
motions . 


5.  The  analysis  predicted  that,  if  the  wing  were  moved  rearward  from 
a  position  of  coincidence  of  wing  quarter  chord  and  missile  center  of 
gravity  (no  horizontal  tail  present),  the  speed  of  flutter  involving  body 
motions  would  first  decrease  sharply  until  a  certain  wing  position  rela¬ 
tive  to  the  center- of-gravity  location  was  reached  and  then  would  increase 
slowly  as  the  wing  was  moved  farther  to  the  rear. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  September  28,  1950. 
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TABLE  II 

MISSILE  PARAMETERS 


[Missile  center  of  gravity  is  17  inches  ahead  of  wing 
leading  edge;  missile  mass  moment  in  pitch  is 
14.8  slug-ft^;  missile  weight  is  236  lb;  air 
density  p  is  0.002272  slug  per  cu  f 1 7j 


Type  of 
analysis 

Degrees  of 
freedom  included 

_ m _ 

Flutter- speed 
coefficient, 

v/bo^i 

Flutter 

frequency  ratio, 

“Mi 

Root  B 

(3) 

Root  A 

Root  B 

a 

F 

T 

P 

V 

9.6 

22.5 

0.26 

3.40 

b 

F 

T 

- 

- 

— 

22.5 

— 

3.40 

c 

F 

- 

P 

- 

7.3 

— 

0.20 

— 

d 

F 

T 

P 

- 

6.4 

22.5 

0.17 

3.40 

e 

F 

- 

P 

V 

11.2 

— 

0.287 

— 

^Designation  of  degrees  of  freedom 
F  -  wing  flexure  (first  bending) 

T  -  wing  torsion  (first  torsion) 

P  -  rigid-body  pitching 
V  -  rigid-body  vertical  translation 

^Root  A  -  critical  body- freedom  root 
3Root  B  -  critical  wing  bending-torsion  root 
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OI25456 


Time  from  launching,  sec 

Figure  6.-  Variations  of  longitudinal  acceleration  with  time. 


Figure  7*-  Portion  of  telemeter  record  of  model  E 


Root  B  (bending  torsion  flutter) 


i  o  ^  s 

M  £  H  K 
§  d  UU 

9  .3  .<< 

O  J  p-  £  S5 

h4  H  Sfc 


co 

cp 


s* 


n4  H  B  fc 


■g 

£ 


< 

u 

•< 

55 


